Interior estimates for solutions of Abreu's 

equation 



S. K. Donaldson 



February 1, 2008 



1 Introduction 

In this paper we study a nonlinear, fourth order, partial differential equation 
for a convex function u on an open set f2 in R™. The equation can be written 

as 



Here (u^ ) denotes the inverse of the Hessian matrix «y = d ®iQ x j ■ We call this 
PDE Abreu's equation since the expression S(u) appears in [1], in the study of 
the differential geometry of toric varieties. In this paper we will be primarily 
interested in the case when A is a constant. Solutions to Equation 1 then 
correspond to certain Kahler metrics of constant scalar curvature, as we will 
recall in more detail in Section 5 below. Our purpose is to derive a priori 
estimates for solutions of Abreu's equation, which can be applied to existence 
questions for such constant scalar curvature metrics, on the lines of [6]. However 
in the present paper we will keep the differential geometry in the background, 
and concentrate on the PDE aspects of the equation. 

Abreu's equation can be fitted, as a limiting case, into a class of equations 
considered by Trudinger and Wang [12], [13]. These authors study the Euler- 
Lagrange equation of the functional 



(or this functional plus lower order terms). The Legendre transform interchanges 
the equations with parameters a and 1 — a. There are two exceptional cases, 



S(u) = A 

where A is some given function and S(u) denotes the expression 




(1) 
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when a = or 1, when the Euler-Lagrange equations are trivial. Abreu's 
equation, which is the Euler-Lagrange equation associated to the functional 



is the natural limit of the Trudinger-Wang family when a — > 0. Indeed, the 
Trudinger-Wang equations (in the absence of lower order terms) can be written 
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We will study Abreu's equation augmented by some specific boundary con- 
ditions. These depend on a measure a on the boundary of fl. We will consider 
two cases 

Case 1 £1 is the interior of a bounded polytope, defined by a finite number of 
linear inequalities with n codimension-1 faces of dQ meeting at each vertex. The 
measure a is a constant muptiple of the Lebesgue measure on each codimension-1 
face of dQ. 

Case 2 fl is a bounded domain with strictly convex smooth boundary and a is 
a smooth positive measure on c?fi. 

The first case is the one which is relevant to toric varieties and is our main 
concern. We include the second case because it seems to lead to a natural PDE 
problem. In either case we define a class of convex functions <Sn i<7 satisfying 
boundary conditions depending on a. The detailed definitions are given in 
Section 2.2 below, but roughly speaking we require that a function u in Sa,a 
behaves like o~ x d log d where d is the distance to the boundary and a is regarded 
as a function, i.e. a multiple of the area measure on 90. (Note here that 
the whole theory is affine- invariant, and does not depend on the choice of a 
Euclidean metric on R n , but for simplicity we will sometimes, as just above, 
express things in terms of the metric structure, although this is not playing 
any real role.) In this paper we study solutions u of Abreu's equation which 
lie in <Sn jCr . As explained in [6] and in Section 2.2 below, these arise when one 
considers the problem of minimising the functional 



over the set of smooth convex functions on Q, with L 1 boundary values. It is 
likely (although this requires proof) that any extrema of this functional lies in 
our space Sn,a- 

As explained in [6], simple examples show that for some data Q,A,da there 
are no solutions of Abreu's equation with the given boundary behaviour. By 




as 






2 



the same token, one needs further assumptions before a priori estimates of a 
solution can be obtained. The condition which we expect to be appropriate 
involves the linear part of the functional T; 

£(/) = / fd° - [ fAdfi. (3) 

Fix a base point in the interior of ft and call a function normalised if it vanishes, 
along with its first partial derivatives, at this base point. We consider the 
following condition on the data (ft, A, da): 

Condition 1 The functional C vanishes on affine-linear functions and there is 
some A > such that 

C{f) > A" 1 f fda, 

JdQ. 

for all convex normalised functions f on ft. 

Of course when this condition holds we can fix A — A(ft, a, A) by taking 
the best possible constant. It is shown in [6], at least in the case when the 
dimension n is 2 and ft is a polygon, that this is a necessary condition for the 
existence of a solution. Our goal in this paper is to derive a priori interior 
estimates for solutions assuming this condition. To streamline the statement we 
introduce the following terminology. We say a function C(ft, a, A, A, d) where A 
and d are positive real variables is tame if it is continuous with respect to the 
natural topology on the space of variables (ft, a, A, A, d) (We use the C^ c H L°° 
topology on A; in Case 1 the space is divided into components labelled by the 
combinatorics of the faces and in Case 2 we use the C°° topology on ft,cr). 
We use the same terminology for functions that depend on some subset of the 
variables. Our main result is 

Theorem 1 Suppose the dimension n is 2. There are tame functions K,C P 
for p = 0, 1, . . . such that if A is a smooth bounded function in ft, and (ft, a, A) 
satisfies Condition 1, then any normalised solution u in <Sq i(7 of Abreu's equation 
satisfies 

K' 1 < (Uy) < K 

and 

W p u\ < C p 

where the argument d is the distance to the boundary of ft and the argument A 
is A(ft, a, A) 

Of course we need not take the definition of "tame" functions C(ft, a, A, A, d) 
too seriously. With a little labour we could make all of our estimates com- 
pletely explicit. The definition is tailored to the continuity method: if we have 
a continuous 1-parameter family of such problems defined by data (Q t ,a t ,A t ) 
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with solutions for t < to then provided A (fit, at, A t ) stays bounded the solutions 
cannot blow up in the interior as t — > to- 

We present a variety of different arguments to establish these interior esti- 
mates. We hope this variety is justified by the desire to extend the results, in the 
future, in various directions: to the behaviour near the boundary and to higher 
dimensions. In Section 2 we bring together some more elementary preliminaries 
and in particular show that Condition 1 gives C 1 bounds on the solution. After 
this, the crucial intermediate goal is to obtain upper and lower bounds on the 
determinant of the Hessian (uij ) . In Section 3 we consider the case when the 
dimension n is 2 and A is a constant. We find a special argument in this case 
using a property of solutions of general elliptic equations in 2 dimensions. In 
Section 4 we find other arguments, using the maximum principle, which apply 
in any dimension. We get a lower bound on the determinant in all cases and, 
using the technique of Trudinger and Wang, an upper bound involving also a 
"modulus of convexity". As we explain in (5.1), in dimension 2, this modulus 
of convexity is controlled by the lower bound on the determinant, using an old 
result of Heinz. Once we have established upper and lower bounds on the de- 
terminant we can appeal to sophisticated analysis of Caffarelli and Gutierrez to 
complete the proof of Theorem 1. This is explained in (5.1). In the case when 
A is a constant and ft is a polygon in R 2 , which is our main interest, we give an 
alternative proof in the remainder of Section 5. This avoids the deep analysis 
of Caffarelli and Gutierrez (and perhaps gives more explicit estimates) but uses 
their geometric results about the "sections" of a convex function in an essential 
way. The other ingredients are L 2 arguments, which make contact with Kahlcr 
geometry, a variant of Pogorelov's Lemma and standard linear theory (Moser 
iteration). We also explain that, in this case, one can avoid using Heinz's re- 
sult, by combing results from Sections 3 and 4. In sum, in the case when A is 
constant and fi C R 2 is a polygon we get a self-contained proof of Theorem 1, 
assuming only Chapter 3 of [9] , and material from the textbook [7] . 

2 Preliminaries 

2.1 Miscellaneous formulae 

This subsection consists of entirely elementary material. We will give a number 
of useful equivalent forms of Abreu's equation, which are obtained by straight - 
foward manipulation. Throughout we use traditional tensor calculus notation, 
with summation convention. 

Begin by considering any convex function u defined on an open set in R". We 
write for the Hessian, u l i for its inverse and U % i for the matrix of cofactors 
i.e. — det(iiij)u l: > . We can associate to u two second order, elliptic, linear 
differential operators 




(4) 
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Q(f) = {U ij fi) j ■ (5) 
Let L = logdet(ujj). Then we have identities 

U=u ah u abl , (6) 

uf = -u^u aM u bk . (7) 

Thus 

«%•=-«?• (8) 

This means, first, that 

Uf = {e- L u ij ). = t 1 (uf - ,/'••/,,) = 0. 
Hence the operator Q can be written as 

Q(f) = u ij f ij . (9) 

Now the first form of Abreu's equation, as in the Introduction is 

4 = ~ A - 

We define a vector field v — (w- 7 ) by 

= -uf (10) 

So our second form of Abreu's equation is 

v] = A. (11) 

(Of course, the left hand side of this expression is the ordinary divergence of 
the vector field v.) On the other hand, by Equation 8, the vector field can also 
be expressed as 

v j — u 10 Lj, 

so we have our third form of Abreu's equation 

(v**Li). = A (12) 

that is 

P(L) = A. 
Expanding out the derivative we have 

(u ij Li). = u ij L ii:j + ujL, = u ij {Lij - LiLj) , 

so we get our fourth form of the equation 

u ij {L ij -L i L j )=A. (13) 



Finally, if we write 

F = det(uy) -1 = er L , 

then 

u^Fij = u i: >{L l0 - LiL^e' 1, = -Ae~ L , 
so we get our fifth form of the equation 

Q(F) = —A. (14) 
2.2 The boundary conditions. 

We begin by giving a precise definition of the set Sq j(7 of convex functions on 
£1. The definitions arc different in the two cases. We start with Case 1, when 
ft is a polytope. For any point P of dfl we can choose affine co-ordinates Xi on 
R™ such that P has co-ordinates Xi = 0, i = 1, . . . n and a neighbourhood of p 
in ft is defined by p inequalities 

Xl,X2, ...,X p >0. 

We can also choose the co-ordinates so that the normal derivative of x,- L on the 
face xi = of the boundary is a -1 . We call such co-ordinates adapted to at 
P. 

Definition 1 In Case 1 the set Sa,a consists of continuous convex functions u 
on ft such that 

• u is smooth and strictly convex in fi, 

• The restriction of u to each face of dQ is smooth and strictly convex; 

• In a neighbourhood of any point P of dil the function u has the form 

U = ^^Xi \0gXi + f 

where X{ are adapted co-ordinates, as above, and f is smooth up to the 
boundary. 

(We say a smooth function is strictly convex if its Hessian is strictly positive.) 

Now turn to Case 2, when f2 has smooth boundary. If P is a point of dil we 
can choose local co-ordinates £, rji, . . . f] n -\ near P so that d£l is given by the 
equation £ = and the normal derivative of £ on the boundary is cr _1 . Again, 
we call such co-ordinates adapted. Define functions a p of a positive real variable 
by 

ai(t) = -logt, a 2 (t) = t- 1 ,a 3 {t)=t- 2 . 

Definition 2 In Case 2 the set <Sq ;(7 consists of continous convex functions u 
on ft such that 
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• u is smooth and strictly convex on f2; 

• In a neighbourhood of any point P of dfl there are adapted co-ordinates 
(£,7?i) in which 

u = £log£ + / 

where for p>l,p + q<3 

|Vj?V«/|=o(a p (0). 

as £-> 0. 

Here the notation V£ 4 V2 means any partial derivative of order p in the 
variable £ and total order q in the r\i. 

Now, for small positive 5, let ilg C be the set of points distance at least 6 
from <9f2, i.e. a "parallel" copy of <9f2. If u e <Sji !(T and \ is a smooth function 
on Q, integration-by-parts over fig gives the fundamental identity: 

/' t^Xij = [ 4* - [ -^X + ^XJ- (15) 
Jn s Jn s Jdn s 

We can write the boundary terms as 

X "norm + VxX, 



J 

Jdn 5 



where v is the vector field introduced in (10) above, w n orm is its normal compo- 
nent and X is the vector field, defined on a neighborhood of <9f2 by 

X j = u ij Vj, 

Vj being the unit normal to d£l at the closest boundary point. (In Case 1, the 
vector field X will be discontinuous near the "corners" of d$l but this will not 
matter.) 

The main result we need is 

Proposition 1 In either case, if u e Sn,a then as 6 — > 0; \X\ = 0(8) and 
fnorm converges uniformly to a. 

Here \X\ refers to the Euclidean length and we interpret « norm as a function 
on <90 in the obvious way, by taking the closest point of dflg. We assume 
Proposition 1 for the moment. Taking the limit as 6 tends to in Equation 15, 
we obtain 

Corollary 1 Let u he in Sq !<7 with S(u) = —uV € L°°(0) and let x oe a 
continuous, convex function on Q, smooth in the interior and with Vx = o(d~ 1 ), 
where d is the distance to dfl. Then u^Xij is integrable in il and 

ulJ Xij = u ijX+ Xd<7- 
Jn Jon 
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The main application wc make of Corollary 1 is the case when \ = u - 
It is clear from the definitions that Vu is 0{— logd), hence o(d^ 1 ), near the 
boundary, and UijU^ = n. So we have the identity 



where A = —ufj and C is the linear functional defined in (3). Thus we obtain 

Corollary 2 Suppose A e L°°({V) and that C satisfies Condition 1. Then if u 
is a normalised function in <Sq i(7 which satisfies Abreu's equation S(u) = A we 
have 



A simple argument ([6] , Lemma 5.2.3) shows that the integral over the boundary 
for a normalised convex function, controls the derivative in the interior and we 
have 

Corollary 3 Under the same hypotheses as Corollary 2, 



where C depends tamely on (ft, a, A). 

This first derivative bound is the seed which we wish to develop in this paper 
to obtain bounds on higher derivatives. 

We mention some other applications of Corollary 1 . Here we will restrict for 
simplicity to the case when f2 is a polytope. 

Corollary 4 In the case when Q C R" is a polytope, for any u G 5n jCr , 
| logdetwjjl is integrable over Q so for any A e L°°(f2) the functional Ta is de- 
fined on <5>n,<r- The functional Ta is convex on Sn,<j and the equation S(u) = A 
has at most one normalised solution u in <So ;Cr . For such u, 



where C is a tame function of SI, a, ||A||i,oo, A. 

The proof of this Corollary follows easily from the results in ([6], subsections 
(3.3) and (5.1)). The restriction to Case 1 arises because in this case, as is will 
be clear from Proposition 5 below, S(u) is bounded for any u e So,<r- In Case 
2 we have not analysed the behaviour of S(u) near the boundary, for general 
elements of Sq i0 -, and we leave this issue to be discussed elsewhere. 



Cu = n 




u da < n\. 



|Vu| < Cd 



>—n 
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2.3 Proof of Proposition 1. 

In the course of the proof we will also establish some other properties of functions 
in <Sn ;Cr . Before beginning it is worth pointing that, in a sense the proofs of 
these boundary properties need not be taken too seriously at this stage of the 
development of the theory. At this stage we are in essence free to choose the 
definition of Sq^ and we could impose any reasonable conditions we choose. 
For example we we could take the conclusions of Proposition 1 as part of the 
definition of Sn,a- The discussion will only acquire an edge when one goes on 
to the existence theory. It is possible that the detailed definitions may need to 
be modified then. At the present stage, the definitions we have concocted serve 
to indicate at least the nature of the solutions we want to consider while being, 
we hope, sufficiently general to permit a sensible existence theory. 

We begin with Case 1. This is not very different from the analysis of metrics 
on toric varities in [1], [8], but we include a discussion for completeness. We 
work in adapted co-ordinates around a boundary point, so 

v 

u = ^2 x l log x l + f 

where / is smooth up to the boundary. By the second condition of Definition 1 
we may choose the co-ordinates so that at x = 

d 2 .f = ^ 
dxidxj 4 "" 

for p + 1 < i, j < n. Thus 

(Uij) =diag(x^ 1 ,...,a;p 1 ,l...,l) + (-0 y -), (16) 
where tpij are smooth up to the boundary and tp vanishes at x = for i, j > p+1. 

Proposition 2 • det(uy ) X ~\_ • • ■ X rp A where A is smooth up to the bound- 
ary and A(0) = 1. 

• There are functions fi, g-ij, hij, all smooth up to the boundary and with 
fi(0) = 1, hij(0) = 0, such that 

{Uij) = diag(/ia;i, . . . , f p x p , f p+1 , . . . f n ) + (V 3 ) , 

where 

for 1 < i, j < p; 
for 1 < i < p, j > p; 
for i,j > p. 



G ^ — X^X j Cfij 



a 11 = hi 
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It is completely straightforward to verify Proposition 1, in Case 1, given this 
Proposition 2. Notice that, according to Proposition 2, the matrix u IJ is smooth 
up to the boundary. 

Let A be the diagonal matrix 

A = diag(zl /2 ,... a; y 2 ,l...,l). 

To prove the first item of Proposition 2 we write H for the Hessian matrix (uij ) 
and consider the matrix A 2 H. This has the form 1 + E where the entries Eij of 
the matrix E are smooth up to the boundary and E^ vanishes at except for 
the range 1 < j '• < p, i > p. Thus, at 0, the matrix E is strictly lower-triangular 
and so A = det(l + E) is a smooth function taking the value 1 at x = 0, and 

det(uy ) = det A -2 A = x^ 1 . ..x^A. 

To prove the second item we consider the symmetric matrix AHA. We write 
A_ffA = 1 + (Fij) and Equation 16 yields 

J? — A V 2 1 / 2 
fij — AijX i X- 



for 1 <i,j<p; 



Fij — Fji — BijxJ , 



for 1 < i < p and j > p where Aij , Bij are smooth up to the boundary. In 
the remaining block, p + 1 < i, j < n, the F^ are smooth up to the boundary, 
vanishing at 0. Let ad(AiJA) be the matrix of co-factors, or "adjugate matrix", 
so that 

(AHA)- 1 = dct(AiJA)" 1 ad(A J ffA). 

We claim that ad(A_ff A) has the form 1 + F^ where F' satisfies the same con- 
ditions as F above; that is, 

1 ' ij ~ A ij X i X j 

for 1 < i,j < p; 

for 1 < i < p and j > p etc. The proof of this is a straightforward matter of 
considering the various terms in the cofactor determinants which we leave to 
the reader. Given this, the second item of the Proposition follows by writing 

(u ij ) = H- 1 = A' 1 A ad(Atf A) A. 

We now turn to Case 2. For simplicity we will consider the case when 
n = 2 (see the remarks at the beginning of this subsection). Moreover, to make 
the calculations easier, we will consider a special kind of adapted co-ordinate. 
Suppose the point P is the origin and let the boundary of SI be represented by 
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a graph xi = Q(x 2 ) where Q(0) = Q'(0) = 0,Q"(0) < 0. Then we take the 
adapted co-ordinates 

V = x 2 , £ = p{f]){Q{x 2 ) - £1), 

for a positive function p equal to cx _1 , where <r is regarded as a function of £2 
using the obvious parametrisation of the boundary. Without loss of generality 
we can compute at points on the x 2 axis, where rj = 0. We write d Xl ,d X2 
for g^j-, and d$, d v for -§^,-§^- Then, transforming between the co-ordinate 
systems (x\,x 2 ) and (£,77), we have 

9 Xl = P®z , d X2 =d v + Fd£, 

where F(£, rj) has the form F = £A(r]) + B(r]) and 5(0) = 0, B'(0) < 0. Thus 

d 2 Xl = P 2 dl (17) 
d Xl d X2 = P 2 d^+pFdl + pAd i7 (18) 
8 2 X2 - (^ + 2F0 £ ft, + F 2 dl) + (£A' + B' + A(£A + B))d s . (19) 

Applying this to a function u = (£ log £ — £) + / in <Sn i<7 (where we have replaced 
/ by / — £ in Definition 2, which obviously makes no difference) we get 

1 ^ ~ V B'log£ y ^ /3 7 
where 

a = ^/, /3 = pAlog£ + 9 Xl 9 X2 /, 7 = (^ , + B' + ^A + B))Ioge + ^ 2 /. 
Then 



where 



Write 



det(uij) = (p 2 B')^~ 1 log£ A, 



pV V s'logey /9 2 s'iog£ 

a£ 7 /?£ 



Then 

and the inverse matrix is 



p 2 ^ B'logC Byiogf 
A = (1 + A)(1+m)-/^, 
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Now first we evaluate at a point where rj = 0, so B = 0. Then we have 

a = p 2 d i d ( f, 

and by the definition of Sn,a this is o(£ _1 ). So A is o(l). We have 

(3 = P 2 d^ n f + rfAd&f + pA{d t f + logO, 

and applying the definition of Sn.a to the various terms we see that (3 = 
0(| log ^|). This means that v is O(0 and [3v is 0(£|log£|), which are both 
o(l). Similarly, applying the definitions, we find that p is o(l). We see from 
this first that the vector field Z, which has components u n 7 u 12 is 0(0, so 
verifying the first item of Proposition 1 . We also see that 

det(««)~! 7 r 1 |k>g£|, (20) 
since A <~ 1. To complete the proof we need to differentiate again. We have 

L = logdet(uy ) = - log£ + log(- logO + log(-B') + log/? 2 + log A, 
Thus, when r) = 0, 

d Xl L = piC 1 + (HlogO- 1 ) + P(^), (21) 
8 X2 L = ^L+^. (22) 

We claim that d v A is O(l) and d^A is o(^ _1 ). Given this, we have 

v 1 = A-^p-^i +p) + (0), « 2 = a-^-^- 1 + od log^r 1 ), 

and we see that the normal component v 1 converges to a = p^ 1 as desired. 
To verify the claim we have to show that d v \,d v p,d v ((3v) are all O(l) and 
d^X,d^p,d^((3v) are all o(^ _1 ). This is just a matter of differentiating the for- 
mulae defining A, p, v and applying the definition of iSq i<7 to each term: we omit 
the details. 

3 The two dimensional case. 
3.1 The conjugate function 

Throughout this Section 3 we suppose that A is a constant and u is a solution 
of Abreu's equation in ft. We begin with second form of the equation 

div(v) = v\ = A. 
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The "radial" vector field x % on R" has divergence n, so if we define a vector 
field w by 

w l = v i - -x\ (23) 
n 

then w\ — 0: the vector field w has divergence zero. Now define a function h by 

h = u-u k x k . (24) 

Then 

hi=Ui- u ki x k - u k 8 k = -u ki x k , (25) 

so 

u ij hi = -x j . (26) 

Thus, using Equation 8, we have 

w l = u ij L h (27) 

where L = L + £h. 

We now specialise (for the rest of this subsection) to the case when n = 2. 
The special feature here is that divergence-free vector fields can be represented 
by Hamiltonians, so there is a function H (unique up to a constant) with 

w l = e ij Hj, 

where e y is the skew-symmetric tensor with e 12 = 1. 
Lemma 1 The function H satisfies the equation Q{H) = 0. 
To see this we write 

L b = -u bk e kj Hj. 

Then we have 

L ab e ab = 

by the symmetry of second derivatives, whereas, differentiating Equation 27, 
L ab e ab = - (e ab Ub k e kj Hj) a = Q(H), 

since U aj = e ab u bk e kj . 

We call H the conjugate function since the relationship between the functions 
H and L, is analogous to that between conjugate harmonic functions in two 
dimensions, except that they satisfy different linear elliptic equations: P(L) = 
0,Q(H) = 0. 

Next we use the boundary conditions, so we suppose that u is in Sa yCr . We 
have shown in Proposition 1 that the normal component of the vector field v 
on the boundary can be identified with given measure a. Likewise, the vector 
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field ^-x 1 defines a (signed) measure dr on the boundary. The condition that C 
vanishes on the constants implies that 



Vol(<9ft, a) = AVol(fi) 



(28) 



which in turn implies 




da — dr = 0. 



Thus there is a function b on dfl, unique up to a constant, with db = a — r. 
Observe that, in the case when f2 is a polygon, the function b is linear on each 
face of dQ. For in this case da is constant multiple of the Lebesgue measure on 
each face and the normal component of the radial field is also constant on each 
face. 

Lemma 2 The function H extends continously to with boundary value b (up 
to the addition of a constant). 

First suppose that H extends smoothly up to the boundary. The relation 
w 3 — e u jj. asserts that the normal component of w is equal to the tangential 
derivative of H . But the normal component of w is the derivative of b, so the 
derivative of H — b vanishes. For the general case we apply the same argument 
to a slightly smaller domain and take a limit, the details are straightforward. 

We can now apply a special result ([7], Lemma 12.6) for solutions of elliptic 
equations 



in two dimensions whose boundary values satisfy a "three point" condition. 
Here we use Equation 9 to see that the equation Q(H) = has this form, 
with a 4 - 7 = U % K The three point condition requires that for any three points 
X\,Xi, X-i on dfl the slope of the plane in R 3 containing (Xj, H{Xi)) is bounded 
by some fixed K . This holds in our situation because the restriction of H is b, 
which is smooth in the strictly convex case and lincar-on-faces in the polygonal 
case. In fact, a little reflection shows that in Case 1 we can take K to be 
the maximum of the slopes attained by taking X\, X2, X3 to be vertices of the 
polygon. So we can apply the result of [7] quoted above, and deduce that the 
derivative of H is bounded by the constant throughout SI. Thus the vector field 
w satisfies an L°° bound over Q and hence the same is true of v. To sum up we 
have 

Theorem 2 Suppose u <G <So ;Cr is a solution of Abreu's equation where A is 
a constant and the dimension n is 2. Then there is a constant K, depending 
tamely on Q 7 a, such that \v\ < K throughout Q. 
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3.2 Integrating bounds on v 

In this subsection we study the implications for L = log det Uij of an L°° bound 
on the vector field v % = —uj associated to a convex function u on a convex 
domain D C R™. We prove two results, under different hypotheses. For the 
first, recall that the derivative Vu is a diffeomorphism from D to its image 
(V«)(D) C R™. 

Lemma 3 Suppose that the image (Vu)(D) is convex and \v\ < K on D. Then 
for any x,y e D 

\L{x)-L{y)\<K\(Vu)(x)-(Vu)(y)\. 

To see this, let & be the standard Euclidean co-ordinates on the image (Vu)(D). 
(More invariantly, Vu should be thought of as mapping D to the dual space 
(R™)* and & are just the dual co-ordinates of x % .) In traditional notation, 

& = u t so 

dxj ~ Uij ■ ,K, " ■ 
Now v % = u^Lj (see (8)) but by the chain rule 

In other words, if we define a function L* on (Vu)(D) to be the composite 
L o (Vu)" 1 then v is the derivative of L* and hence |Vi*| < K. If (Vtt)(£>) 
is convex then Vu(x), Vw(j/) can be joined by a line segment in Vw(-D) and 
the result follows immediately by integrating the derivative bound along the 
segment. 

The next results applies to solutions of Abreu's equation, but avoids the 
convexity hypothesis on the image under Vu. We write Avjy(L) for the mean 
value 

Av D {L) = (Vol(D))- 1 f L. 

J D 

Theorem 3 Suppose that D C R" is a bounded convex domain with smooth 
boundary and xq G D is a base point. Let R be the maximal distance from Xq to 
a point of dD. Let u be a solution of Abreu's equation S(u) = A in D, smooth 
up to the boundary and normalised at x . Then 

\L(x ) - Av D (L)\ <C I udv, 

JdD 

where dv is the standard Riemannian volume form on dD and 
C = (nVoliD))- 1 (nsupH +i?sup|A 

V D D 
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To prove this, we can suppose that xo is the origin and let (r, 9) denote standard 
"polar-coordinates" on R™ \ {0} (so 9 £ 5 rl_1 ). Suppose that dD is given in 
these co-ordinates by an equation r = p(9_). Let / be the function 

/-('-^) 

on R™ \ {0} and define a vector field ( on R™ \ {0} by 

C = M 

Then, away from the origin, 

o n d p n o n d p n 

C = x l (l - —) - x 1 ——— = nil - —) - r—— = n 

and C, vanishes on the boundary of D. It is clear that £ satisfies the distributional 
equation 

Q ^n-nVol(D) So, 

where So is the delta-function at the origin. Applying this distribution to the 
function L we get 

n(Vol(£>)L(0) - / L)= f QU. 

J D J D 

Now recall that h = u — UiX 1 satisfies hi = —UijX J . Then Li = UijV J , and 
x l Li — UijX l v^ = —hjV J = — (hv 3 ). + hvj. 

Using Abreu's equation in the form v 3 - = A we have 

x l Li = - (hv j ). + Ah. 

Now 

/ CU = f fx' Li = - I f{hiP)j + f Afh. 

JD JD JD JD 

We can integrate by parts on the first term to get 

/ CU= f (hf jV i+Afh). 

JD JD 

There is no boundary term from dfl since / vanishes on the boundary. There is 
also no extra term caused by the singularity of / at the origin since h vanishes 
to second order there. To sum up, we have the identity, 

nVol(D) (£(0) - Av D (L)) = {I) + (II), (29) 
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where 



(I) = / hf s v> (II) = / Afh. 

J D J D 

We next estimate the size of the integrals (7) and (77). We have 



df _ p n 1 df _ p™" 1 dp 

so 

\\7f\ = - 

r n+l 

where 



iv/l = 



Now let 7f = sup D \v\ so we have 

\fy\<^^w(0_). 



Thus 



|(/)| < nK J D W(9)p(er-i-L- = 2K J D P n - l W(6)^drd6_, 

in an obvious notation. Now consider this integral along a fixed ray 8 
constant. In polar co-ordinates we can write 

du 

h=u-r—, 
or 



so 



. I . _o — 1 — 9 j 

— (r u) = —r u + r — = —r h. 
or or 



So 

ff(£) 



J ^dr = p- 1 u(p(6),e)-e- 1 u(e,6), 

and the term from the lower limit tends to zero with e since the derivative of u 
vanishes at the origin. So we have 

\{I)\<nK [ P (&) n - 2 W{6)ud6. (30) 
JdD 

For the other integral we have —ru r < h < 0, and / < so 

<(//)< sup \A\ f ^■ru r r n - 1 drd9 = sup \A\ [ p(9) n ud9. (31) 

J D r ™ JdD 
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We now transform the integrals over dD to the Ricmannian area form dv. 
By straightforward calculus 



dv = p 11 - 1 ^ 1 + \^-=j d9 = p n - 2 W(8)d9. 
So inequality 30 is just 




\{I)\<nK [ udv. 



JdD 

Also, dv > p n ~ 1 d9 so inequality 31 gives 

|(7/)| < sup \A\ I updv < Rsup \A\ / udv. 



JdD JdD 



and putting these together gives the result stated. 

We now apply these results in the case when n = 2 and A is a constant, 
using Theorem 2. We obtain 

Theorem 4 Let CI C R 2 and u e Sn.a be a normalised solution of Abreu's 
equation with A constant. Suppose Condition 1 holds. Then 



where d is the distance to the boundary of CI and Co , C\ , Ci depends tamely on 
CI, a, A. 

We can prove this using the simpler result Lemma 3. Any function u in 
is continuous on CI and so attains a minimum value on O, but it is clear from 
the definition of <Sq, ct that this cannot be attained on the boundary. Changing 
u by the addition of any linear function we see that for u G Sa,a the image 
(Vw)(f2) is the whole of R™. Let B be a small disc about the base point. We 
know that Vu is bounded on B by Corollary 3, and this gives a bound on the 
area of Vu(i?) , which is 



So there is a c, depending tamely on the data, such that we can find some point 
in B where L < c. Now the result follows by combining Corollary 3, Theorem 
2 and Lemma 3. 

We can give a slightly different proof of the above result using Theorem 3 in 
place of Lemma 3. In the case when CI C R 2 is a polygon and A is a constant 
we can also obtain an interior lower bound on L by applying Corollary 4, which 
gives a bound on 



However we will not discuss this in detail, since we will get a better lower bound 
in the next section. 



L<C + C!\Vu\ <C 2 d- 2 , 
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4 Applications of the maximum principle 

4.1 Lower bound for the determinant 

In this Section we will derive a lower bound for L = logdet(uy ), valid in any 
dimension and for any bounded function A. 

Theorem 5 Suppose that u € <Sn, CT is a solution of Abreu's equation for some 
bounded function A. For any a e (0, 1) there is a constant C a (depending only 
on n and a) such that 

detUij > C a d- a (supA) n Diam(ft) 2ll+Q 

throughout il. 

Here we recall that d is the distance to the boundary of £1. Diam (fi) is the 
diameter of fi. 

We can compare this result with the asymptotic behaviour established in 
Section 2 (Proposition 2 and (20)). In Case 1 

detu i3 ;~ dd' 1 , (32) 

as x tends to a generic boundary point (i.e on an (n — l)-dimensional face of 
the boundary) and in Case 2, 

det Uij ~ C 2 eT 1 |logd|. (33) 

To prove the theorem we suppose that we have a function tp on CI which 
satisfies the following conditions. 

1. For each point of fl the matrix 

M i:j = ipij - ^ipj 

is positive definite. 

2. L — ip tends to +oo on dfl. 

Certainly such functions exist. For example, we can take ipe( x ) = e \ x \ 2 f° r small 
e: the second condition holds since I^ooon dil. Given such a function ip, we 
write D = det (M i:7 -). 

By the second condition there is a point p in fl where L — ip attains its 
minimum value. At this point we have 

L t = ip t (Lij - ipij) > 0. 

We take Abreu's equation in the form (Equation 13) 

u ij (Lij - LiLj) = A. 
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At the point p we have 

u ij Lij = A + u ij LiLj = A + u ij 4>iipj, 

and 

u ij Lij > u l3 ipij. 

Thus, at the point p, 

i^tuij - ViVj) = H'-'A/ij < A. 

So u 13 Mij < A, where A = sup fi A. We now use the standard inequality for 
positive definite matrices, 

(det(u tf ) det(M ij )) 1/ " < n^u^M^. (34) 

(This is just the arithmetic-geometric mean inequality for the relative eigenval- 
ues of (uij), (Mij).) Now use the fact that det(u y ') = e~ L to get, at the point 
P, 

e~ L D < (A/n) n , 

or 

L(p) > log D(p)-c, 

where c = n\og(A/n). Then for any other point, q, in Q we have (L — ip)(q) > 
(L - tp)(p) so 

L(q) > L(p) + i){q) - ^(p) > log D(p) + yj(q) - ^(p) - c, 
or in other words 

£(<?)> V(<?) + CV, (35) 

where 

= inf(log_D — ip) — c. 

(Of course, at this stage, could be — oo, in which case Equation 35 is vacu- 
ous.) 

By taking the function ip e , say, we immediately get a lower bound L > const, 
over fl. (This is all we need for our main application below.) To prove Theorem 
5 we need to make a more careful choice of the comparison function ip. In fact 
the optimal choice of this comparison function leads to an interesting Monge- 
Ampere differential inequality which we will digress to explain. We consider, 
for a fixed point q the set of functions tp satisfying conditions (1) and (2) and 
with > — oo. The optimal bound we can get from the argument above is 
given by the supremum over this set of functions -0 of ip(q) + C^. Changing ip 
by the addition of a constant, we may suppose that = 0. In other words we 
have L(q) > X(q) — c, where 

X(q) = sup ip(q) 
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and X is the set of functions which satisfy (1),(2) and in addition 

logD<V- (36) 
This becomes more familiar if we write R = — exp(— tp). Then 

R tl = e-^(?% - ViV'j) = e-^Mij. 

So condition (1) is simply requiring that R be a convex function. We have 
D = e"^ det(-Ry) so Equation 36 becomes 

det(fly) > (-R) n -\ 

We can sum up the discussion in the following way. Given a domain fi, let 
TZn denote the set of negative convex functions R on f2 with 

hm — — = -oo 

x^dn d(x) 

and with 

det(iZy) > (-R) n ~ 1 . 

Define a function pq by 

Pn(x) = - sup R(x). 

ReTln 

Then we have 

Proposition 3 In either Case 1 or Case 2, a solution u e Sq !<7 of Abreu's 
equation with boundary data satisfies 

i , . -i / sup f2 A\ n 
det Uij > p n I — - — I . 

in fl. 

This follows from the argument above and the asymptotic behaviour 32, 
33 (In Case 2 we could strengthen the statement a little, taking account of the 
logd term in 33.) 

We return from this digression to complete the proof of Theorem 5. For this 
we take co-ordinates (xi, . . . , x n -i, y) on R" and consider the function 

on a cylinder Z = {(x,y) : \x\ < 1,0 < y < 1}. Here a <G (0,1) is given, as in 
the statement of the Theorem, and b > will be specified later. We will choose 
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b with b < 2 so that r is negative on Z. Straightforward calculation shows that 

2 ' 

b 



r is convex on Z provided (1 — a)(l — |)> afr. This is equivalent to 



> a, 



2 + 6 

so, whatever a is given, we can choose b so small that the condition holds. Then 



(_ r )n-l 



ay""" (1 - a) 




which is bounded below on Z so det(ry ) > C(— r)™ -1 throughout Z for some 
fixed positive C. Writing R = C~ x r we get a function with det Rij > (— 
on Z and with — R = 0(y a ) as y — > 0. 

Now one checks first that the inequality to be proved is stable under rescaling 
of the domain ft. Then given a point p in our domain f2 we can obviously suppose 
without loss of generality that O lies in the set Z (for some suitable K depending 
on f2) , that the origin lies in dSl and that the minimum distance to the boundary 
of SI from p is achieved at the origin. Then the function R constructed above 
(or, more precisely its restriction to Q) lies in the set TZn, since R/d tends to 
— oo on dfl, and the result follows from Proposition 5. 

It should be possible to sharpen this bound in various ways, and this is 
related to the Monge- Ampere Dirichlet problem for convex, negative, functions 
R on SI : 

det Rij = (-i?)™- 1 R\ dn = 0. 



4.2 An upper bound on the determinant 

In this subsection we will modify the method of Trudinger and Wang in [12], [13] 
to obtain an upper bound on the determinant of (?%■). The argument applies 
in any dimension but the result we obtain requires additional information on 
"modulus of convexity" of u for its application, see (5.1) below. 

Consider a bounded domain with smooth boundary D in R™ and a smooth 
convex function u on D satisfying S(u) = A, with u < in D and w^Oon dD. 
We suppose that u is smooth up to the boundary of D. We consider Euclidean 
metrics gij on R" with determinant 1 and for each such metric let 

C g = m&xg^Uiiij. 

Now let 

C = min C g . 

9 

This defines an invariant C of the function u on D. Another way of expressing 
the definition, is that w n C" is least volume of an ellipsoid containing the image 
of D under the map V« : D — ► R™, where u) n is the volume of the unit ball in 
R™. 
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Theorem 6 In this situation 



in D, where 



(detK)) 1 /" < Q + a_M^ eC ( _ tt) _ s 



a = max(0, — min^4) , M = max(— it). 

D D 



To prove this we let g be the metric with C g — C and consider the function 

/ = — L — nlog(-u) — ag lJ UiUj 

on U, where a is a positive constant to be fixed later. Thus / tends to +00 on 
3D and there is a point p in D where / attains its minimum. At the minumum 
fi = and u^fij > 0. The first of these gives 

777/ ' 

L i = -—^-2ag P9 u p u qi . (37) 

The second gives 

9 

- - Tt ft 

< -u^Ltj - — + -^u tJ UiUj - 2a(u l: > g pq u pi u qj + g pq u p u qlJ ), (38) 

where we have used the fact that u^Uij = n. The crucial step now is to observe 
that 

J gP<l UpjUqj — g Pq U pqi 

the ordinary Euclidean Laplacian of u, in the metric g. We next use the form 
Equation 13 of Abreu's equation to see that 

-/'■'/.,, • '/'■'/-,/-., • A, (39) 

where A = minu S(u). Now Equation 37 gives, at the point p, 

• ■ ( 77-^ o GtTt \ 

u lJ LiLj = u lJ I -^UiUj + 4a g pg g rs u p u r u qi u s:j + 4—g pq UiU p u qj J , 
which can be written in the simpler form 

u lj LiLj = —u^UiUj + 4a 2 g pq g rs u pr u q u s + 4^-^g pq u p u q . 

Our inequalities 38, 39 give 

np 1 ■ - o A:Cxn 

< - A ~(u lJ UiUj) - 4a (g pq u rs u p u q ) g pq u p u q h 

u u u 

+ —{u^u l u )-2ag pq u pg -2au l 3g pq u ptl u q . (40) 
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Now recall that L q = vP k u q jk, so Equation 37 gives 

- ■ ft 

g pq u %: >u plj u q = g pq L p u q = --g pq u p u q - 2au pq u r u s g pr g qs 

This means that three of the terms in (40) cancel. There are also two terms 
involving the expression u^UiUj which we can combine to get 

< —A - — + — fi- (u ij u iUj ) ~ 2—g pq u p u q - 2ag pq u pq . (41) 

Now use the fact that n — n 2 < and 

g pq u p u q < C, 
by the definition of C. Re-arranging, we obtain 

9 pq u pq ^ ( n t ^ f-l\ A 



n \2a I \ u J 2an' 
where A = min A. By the definition of a and M we have 
g^u pq< ( z l\ /n ^ + aM\ 



< — —+C; 

n \ u ) \2a 2a J 

Now, since g pq has determinant 1 we have (as in Equation 34) 

e L = det(utf) < i^^ 1 
so, at the point p, 

L < —n log(— u) + n log P, 
where P = "^° M + C. Going back to the definition of / we obtain 

f(p) > —ctg pq u p u q — nlogP > —aC — nlogP. 

So at any other point of D, we also have/ > —aC — nlogP which gives 

det(uij) 1/n < kC(-u)- 1 , 

with k = exp(^ + logP). Optimising the choice of a one finds that the least 
value of k is 



C\\jm+ — + ^ + l exp \\jm+ — -- , 



where m = \ + Now n min < Ce(m + 2) since ^/m+^<y + l, and this 
gives the inequality stated in the Theorem {i.e. the statement is not quite the 
optimal result: notice that if A > 0, so a — 0, one gets K m i n = 2Ce 1//2 .) 
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5 Estimates for higher derivatives 



5.1 The modulus of convexity and the proof of Theorem 
1 

Suppose we have upper and lower bounds on det(ujj) in the interior of for a 
solution u of Abreu's equation, for any smooth function A and in any dimension. 
We can then obtain estimates on all derivatives, following the discussion of 
Trudingcr and Wang [12], [13], in terms of a "modulus of convexity" of u. This 
uses the results of Caffarelli [3] and Caffarelli and Gutierrez [4] and the method 
exploits the form Equation 14, Q(F) = —A, of Abreu's equation. 

We introduce some notation which we will use more extensively in (5.4) 
below. For a smooth strictly convex function u on any open set U C R™ and 
a point x £ U we let X x be the affine linear function defining the supporting 
hypcrplanc of u at x {i.e. u — \ x vanishes to first order at x). Then we define 
the function H x on U by 

H x (y) = u(y) - X x (y). 

(So H x is the normalisation of u at x, in our previous terminology.) Thus 
H x {y) > with equality if and only if x — y. We can think of the function 
H x {y) of two variables x, y as a kind of "distance function" on U, although it 
need not satisfy the axioms of a metric. For a subset S C U we put H X (S) = 
mf yeS H x (y). 

Now return to the case where u is a convex function on fl and K CC K + 
are compact convex subsets of il with < A < det «y < A on K + . Let 

H(K,K+) = min H x {dK+). 

Caffarelli and Gutierrez prove that there is an a e (0,1) such that for any 
solution / of the linear equation Q(f) = — A and x, y in K 

\f(y)-f(x)\<CH(K,K + y a \x-y\ a 

where C depends only on A, the supremum of |Vit| over K + and the upper and 
lower bounds A, A of det(uy ) over K . [The results of [4] are stated for the case 
when A = but, according to Trudinger and Wang ([13], discussion following 
Lemma 2.4), the arguments go over to the inhomogencous equation.] Thus in 
our situation, taking / = F — det (ibij •) ~ 1 , if we have a positive lower bound 
on H(K,K + ) we get a Holder estimate on F and hence on det(uy ). Then the 
results of Caffarelli in [3] give C 2 ' a bounds on u over K, depending again on the 
H(K, K + ). This gives C a control of the co-efficients of the linearised operator Q 
and we can apply the Schauder estimates ([7] Theorem 6.2) to get C 2 ' a control 
of F, and so on. 

Further, if we have a lower bound on the "modulus of convexity" H(K, K + ) 
we can apply Theorem 6, using the interior bound Corollary 3 on the first 
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derivatives of u, to obtain an upper bound on det(ziy •) over K. Given a point x 
in K we set 

u = H x - l -H{K,K+) 

and 

D = {y : u(i/) < 0}, 

soDc and we can apply Theorem 6 to u to get an upper bound on det(uy ) 
at x. Combining this with the lower bound from (4.1) (and replacing K + by 
K) we can then feed into the the preceding argument. Thus the obstacle to 
proving a result like Theorem 1 in general dimensions n is the need to control 
the H(K, K + ), for interior subsets K CC K + C Q. In dimension 2 we can 
again argue in parallel with Trudinger and Wang in [13]. The result of (4.1) 
gives a lower bound det(wjj) > > 0, over K + , for a solution satisfying our 
boundary conditions. Then a result of Heinz [10], implies that in two dimensions, 

H(K,K + )>C\)£, 

where the constant C depends on K + , the distance from K to the boundary 
of K + and sup K+ |Vu|. Putting all these facts from the literature together, we 
arrive at a proof of our main Theorem 1 . 

As we explained in the Introduction, we will give below an alternative proof, 
in the case when A is constant and Q is a polytope, which avoids the sophisti- 
cated analysis of the Caffarelli-Gutierrez theory but employs some of their basic 
tools. The modulus of convexity will also enter in a crucial way in our argument. 
We close this subsection with two further remarks 

• In the case when A is a constant and is a polygon in R 2 we can combine 
the lower bound Theorem 5 and the upper bound Theorem 4, to obtain 

|Vu| > Cd- a 

near the boundary of f2, for some C, a > 0. This gives easily a lower bound 
on H(K,K + ) for suitable K,K + and means we can avoid appealing to 
the result of Heinz in this case. 

• In two dimensions the main result of Heinz in [10] gives a C 1 '' 3 bound on 
u, in terms of upper and lower bounds for det(uy ). In the case when A 
is a constant we can use Lemma 3, and the bound on the vector field v 
(Theorem 2) to obtain a bound on det(uy ). This gives an alternative 
path, avoiding appeal to [4], but feeding into [3], in this case. 

The strategy for our alternative proof is as follows. In (5.2) we introduce 
two tensors F, G depending on the 4th. order derivatives of the function u and 
related to the Riemann curvature and Ricci tensors of a certain Ricmannian 
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metric. Wc show that our boundary conditions fix the natural L 2 norms of 
these tensors. Then we show in (5.2) that under suitable conditions the L°° 
norm of G (or F) controls the second derivatives of the function u, while in 
(5.3) we show that under suitable hypotheses, including control of the second 
derivatives of u, the L°° norm of G is controlled by the L 2 norm of F. Putting 
together these three ingredients we complete the proof in (5.4), using a scaling 
argument and some of the basic geometrical results of Caffarclli and Gutierrez 
on the sections of a convex function. 

5.2 Curvature identities and L 2 bounds 

For a convex function u on an open set U in R™ we define a 4-index tensor by 

rpab ab 

We can raise and lower indices in the usual way, using the metric mj, setting 

-u u t kl ,h ijk i -u ia Uj b t kl . 
Lemma 4 The tensor F^kl is symmetric in pairs of indices 

Fijkl — FkUj . 

In fact, calculation gives, 

Fijkl = Uijkl - U Xq (ukjqUux + Uik\U 3 l q ) 

which makes the symmetry apparent. 

We now introduce a Ricmannian metric, due to Guillemin [8], on the 2n- 
dimensional manifold U x R™, with co-ordinates r\i in the second factor; 

g = Uijdx l dx^ + dr]idr]j . (42) 

This is in fact a Kahler metric: complex co-ordinates and a Kahler potential 
are furnished by the Legendre transfrom construction. If £j = m are the usual 
transformed co-ordinates we set Zi — & + \/— T f]i so 

dzi = Uijdx-i + \f— 1 drji. (43) 

Lemma 5 The curvature tensor of g is 

-~F ijkl dzidz k ® dzjdzi. 

Of course we can use (43) to express the curvature tensor entirely in terms 
of products of the dx % and drjj , avoiding the Legendre transform. 

The proof of the Lemma is mainly a matter of notation. We use the stan- 
dard fact that if the Kahler metric is expressed by a (Hermitian) matrix-valued 
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function H, in complex co-ordinates, then the curvature, viewed as a matrix of 
2-forms, is —d(H(d(H^ 1 ))). In our case the matrix H has entries 

so r = —H(dH^ 1 ) is the matrix of 1-forms with (A/x) entry 

Now use the fact that = u al -J~ to write this as 

r£ = -u Xv u ai u^idz a , 

and this is just 



So the curvature, as a matrix of 2-forms, has A, fj, entry 



d 

— (ul a ) dzfidz a . 

Expressing the derivatives in terms of the x % variables again, this just means 
that the (1,3) curvature tensor is 

n UZy, 

and lowering an index gives the stated formula for the (0, 4) curvature tensor. 

Of course we can also obtain Lemma 3 via Lemma 4 and the usual symme- 
tries of the curvature tensor, but we preferred to give the direct calculation. In 
fact in what follows we will make little explicit use of the Riemannian metric g, 
and derive most of the formulae we need directly. 

Now define a 2-tensor G by contracting F: 

G\ = F%, (44) 

and likewise G tk ,Gik- If follows from Lemma 3 that the latter are symmetric 
tensors, and contracting in Lemma 4 shows that G is essentially equivalent to 
the Ricci tensor of the metric g. In terms of our vector field v, 

G\ = v k- 

A further contraction yields the scalar invariant S — G\ of u which is of course 
the term appearing in Abreu's equation and which corresponds to the scalar 
curvature of g. 
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It is well-known in Kahlcr geometry that on a compact Kahler manifold the 
L 2 -norms of the scalar curvature, the Ricci curvature and the full curvature 
tensor give essentially equivalent functionals on the metrics in a given Kahler 
class: they are related by topological invariants of the data [5]. In our setting we 
do not necessarily have a compact Kahler manifold available so we will develop 
the corresponding theory directly. 

On account of the symmetry of Lemma 3, the standard square- norm of the 
tensor F, using the metric is 

\Ff = F%FM, 

and this is the same (up to a numerical factor) as the standard square-norm of 
the curvature tensor of g. Similarly 

\G\ 2 = G)G\, 

is essentially the square-norm of the Ricci tensor. We consider a 1-parameter 
family of functions u(t) with 

d 

dt U ^ = £ ' 

so e is a function on U. We write E % i for the ^derivative of w y at t = 0, so 

& = -u ia e ab u bj . 
Proposition 4 The time derivatives at t = satisfy: 

d 



and 



where 



dt (|F| 2 - |G| 2 ) = 2Zj 

l(\ G f-s>) = 2w; 



W l = -E{ k G\ + SE J /. 



J 



The proofs are straightforward calculations. 

Now return to the setting of our convex domain f2 in R™, where we suppose 
we are in Case 1 with a polytope. Any pair uo,«i of functions in Sq can 
be joined by a smooth path u t (for example a linear path) in Sn For 5 > let 

be an interior domain with boundary a distance S from the boundary of fi. 
Then the time derivative of 



2^ 



\F\ 2 
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is the boundary integral 

Z l + W\ (45) 



IdQs 

Lemma 6 The boundary integral (45) tends to with S, uniformly over the 
parameter t G [0,1] 

Given this we obtain 

Corollary 5 In the case when £1 is a polytope there is an invariant x(n,<r) 
such that for any u G Sn,a 



[ \F\ 2 -S 2 =x(n,a). 



J 



It is easy to see that x(f2, a) is a tame function of 0, a. Now if u G Sn t(J has 
S'(u) constant then the value of the constant is fixed by Equation ?? and we see 
that 

J ' Xl ' ) + Vol(O) ' 
is a tame function of tt, a. 

To prove Lemma 5 we go back to Proposition 2. Note first that, in an adapted 
co-ordinate system, the matrix is smooth up to the boundary so the same is 
true for F^j and for the time derivative E*i . (Actually, in Proposition 2 we chose 
a special adapted co-ordinate sytem to diagonalise the Hessian of the function 
in the variables Xi for i > p but it is easy to see that the same conclusions 
hold without this restriction.) Thus the vector field Z, W are smooth up to 
the boundary and we simply need to show that the normal components vanish 
on the boundary. Thus we can restrict to a neighbourhood of a point in an 
(n — l)-dimcnsional face of the boundary, i.e. with p = 1. Now, according to 
Proposition 2, the entries w lj are all products of x 1 with smooth functions so 
u 1 ^ = for all j and k > 1. Thus 

El j =0 k > 1. (46) 

The x 1 derivative of it 11 is fixed by the boundary measure so 

El 1 = 0. (47) 

As for the second derivatives we have 

= =0 fOT M> 1 (48) 

and 

Fji = -v-ji - (49) 
since the boundary measure is constant. It is now completely straightforward 
to check that (46), (47), (48), (49) imply that all the terms vanish in the sums 

Z 1 =-Ei l Ftf +E 1 k 1 G k 1 , W 1 ^-Ej k Gl + SEj\ 
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5.3 Estimates within a section 

We begin by considering a convex function u on a bounded open set, with 
smooth boundary, D C R™. We assume that u is smooth up to the boundary, 
that u < in D and that u vanishes on 3D. Suppose that, at each point of D, 
we have 

u > -c , \G\ < ci, |u| Euc < c 2 , |Vw| Eu c < c 3 , 

for some Co, c\, C2, C3 > 0. Here G denotes the tensor Gij introduced above and 
\G\ is the natural norm computed in the metric defined by u, that is 



On the other hand the quantities |u|euc, |Vw|euc are the norms of the vector 
field v l and the derivative Ui computed with respect to the standard Euclidean 
metric on R™. 

Proposition 5 There is a constant K, depending only on n, Co, C\, C2, C3 such 
that (uij) < _ftT|u| _1 on D. 

The proof is a straightforward variant of Pogorelov's estimate for solutions of 
the Monge- Ampere equation det(wjj) = 1, see [9] Chapter 4. (It is also similar to 
the proof in (4.2) above.) It obviously suffices to estimate the second derivative 
un along the co-ordinate axis, and we set 



Now consider the function log(— u) + f which tends to —00 on the boundary so 
has an interior maximum. At this maximum point we have 



|G| 2 = G)G\ = G ij G ab u ia u jb . 



/ = Iog(«n) + ^. 




(50) 



and 



P(-log(«) + /)<0 



where P{4>) = (ti I3 ^)j. Now 




which gives 



Plog(wn) = 



U J Ullij 



+ 



UjUm 



(51) 




Now 
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Expanding out the derivatives we get 

u lJ Uuij = -ufjUn + u la u jb u lab u li: j - ujuui- (52) 

Here we recognise the expression —u^un as the co-efficient G\\ of the tensor 
G. (fn fact the manipulation above is essentially the familiar identity, in Rie- 
mannian geometry, for A|d0| 2 involving the Ricci tensor, where is a harmonic 
function, but we have preferred to do the calculation directly.) Substituting 
(52) into (51), two terms cancel and we get 

G 1 1 

Plog(un) = 1 u la u jb u lab u li: j y-u^umunj. (53) 

Mil Mil 



Now simple calculations give 



P(Y ) = "ii, (54) 

n-v'u, u ij UiUj 
P(log(-u)) = T J -. (55) 



u 

So we conclude that, at the maximum point, 



+ u n + — + (A) - (B) - (C) < 0, 

u U\\ 

where (A), (B) and (C) are the positive quantities 

(A) = —u la u^u lab u Ujl (56) 

(B) = —u lJ umuuj, (57) 

(C) = ^u ij u iUj . (58) 

We now use the condition (50) on the first derivatives. We may suppose that 
Uij is diagonal at the given maximum point, with diagonal entries un = Xi say. 
Then we have for i ^ 1 

Ui_ _ "Hi 
u Un 

The term (C) is 



We write this as 



(C) = + (C)', 



li 
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where 

i>\ i>l 1 

On the other hand, computing at this point, the other sums become 



( B ) = E^n.' ( 6 °) 

Oncc sees from this that 

(A) (B) - (C)' > 
so we conclude that at the maximum point 

n - u t v l Gn u\ 

1- un H T J — < 0. (61) 

u Uu U Z Un 

This should be compared with the more standard calculation, for solutions of 
the Monge- Ampere equation det(uy ) = 1, when v and G vanish, so two terms 
(61) are absent. 

Now, computing the norm in the diagonal basis, 



so 

I — < a- 

Clearly n — UiV 1 < n + C2C3. Then (61) gives 

— - < 0, 

where N = n + C2C3 + c\Cq. Thus we can adapt the argument from the standard 
case, replacing n by N. If ft. is the function 

u 2 

ft = exp(log(-u) + /) = \u\uu exp(y ), 
we see that at the maximum point for ft 



u 2 

ft 2 -N exp(— -)ft ~ u\ exp(u 2 ) < 0, 
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which gives 
This gives 
where 



h max < exp(^ ( N + \/N* + Au\ 



K 

tin < t r 



^=^exp(|) ^ + v /iV2 + 4 C |). 



5.4 Yang-Mills estimate 

In this subsection we consider a solution of the equation S(u) — A with A 
a constant. Roughly speaking, we show that when the dimension n is 2 and 
given uniform bounds on the Hessian of u, the L 2 norm of the tensor F (the 
"Yang- Mills functional") controls the L°° norm. 

Proposition 6 Suppose u is a convex function on the unit disc D in R 2 , 
smooth up to the boundary, which satisfies the equation S(u) = A where A 
is constant. Set 

£=[ \f\ 2 - 

J D 

If the Hessian of u is bounded above and below 

K- 1 < (uy) < K 



then 

and for any p > 1 



|G(0)| 2 < k(£+£ 3 ) 



[ \F\p < Kp (£ + £ 3 )p 

where k depends only on K and sup D \v\ and k p depends on K and p. 
The proof makes use of the Sobolcv inequalities. In dimension 2 we have 

< C p \\Vcj>\\ L , 

for any p and compactly supported <f> on D (here all norms are the standard 
Euclidean ones). The proof can be modified to give a similar result in dimension 
3 and extended to give information when n = 4 provided £ is sufficiently small, 
in the manner of Uhlcnbeck[14] and, still more, Anderson [2] and Tian and 
Viaclovsky[ll]. In the proof we make more use of the Riemannian metric g on 
D x R 2 defined by the convex function u. 
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The idea of the proof is to exploit the fact that the tensors F and G satisfy 
quasi- linear elliptic equations. To derive these we use the interpretation of 
these tensors as (essentially) the Riemannian curvature and Ricci tensor of the 
metric g see also [11]. (Of course there is no particular difficulty in deriving 
these equations directly, without explicit reference to Riemannian and Kahler 
geometry, but the derivation involves manipulating sixth order derivatives of 
the function u.) 

Lemma 7 Suppose that g is a Kahler metric of constant scalar curvature. Then 
the Riemann tensor Riem and Ricci tensor Ric of g satisfy 

V*VRic = Riem * Ric, 

V* VRiem = Riem * Riem + 2V'V"Ric, 

where * denotes appropriate natural algebraic bilinear forms and V', V" are the 
(1,0) and (0,1) components of the covariant derivative. 

To derive these identities we can consider more generally the curvature tensor 
$ of a holomorphic vector bundle E over a Kahler manifold M. Then we have 
9-operators 

8 : (EndE) -» W< q+1 (End£) . 

The Laplacians V* V and = 2(9 d + dd*) differ by a Weitzenbock formula 
involving the curvature of the bundle and base manifold. Since 9$ = 0we have 

= 2dd*<S>. 

Now the Kahler identities give 

9*$ == iV'(A$), 

where A : il 1 ' 1 (EndE) — ► O (End_E) is the trace on the form component. Thus 
V*V$ is equal to 2«<9V(A$) = 2iV"V'(A$) plus a bilinear algebraic term 
involving <f> and the curvature of the base manifold. 

We apply this first to the bundle E = ATM, the anticanonical line bundle. 
In this case $ is essentially the Ricci tensor and A$ is the scalar curvature. 
So the constant scalar curvature condition gives Ag<I> = and the first formula 
follows from the discussion above. The second formula does not use the constant 
scalar curvature condition. It follows from the discussion above taking E = TM , 
when $ is the Riemann curvature tensor and A$ is the Ricci tensor. 

In our situation we deduce that F, G satisfy equations which we write, rather 
schematically as 

V*VG = F*G (62) 
V*VF = F*F + 2V'V"G (63) 
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The strategy now is to take the L 2 inner product of these identities with 
suitable compactly-supported tensor fields. It is convenient here to regard g as a 
metric on the 4-manifold D x T 2 where T 2 is the torus R 2 /Z 2 . However all our 
data will depend only on the D variables and the torus factor plays an entirely 
passive role: the push-forward of the volume form on the 4-manifold to D is 
the standard Lebesgue meausure and the Laplacian of g, applied to T-invariant 
functions is the operator P. Notice the crucial fact that for fixed K the metric 
g is uniformly equivalent to the standard Euclidean metric. 

To prove the first part of Proposition 6 we fix a standard cut-off function 
on the disc, equal to 1 in a neighbourhood of 0, and take the L 2 inner product 
with 2 G on either side of (62). This yields 

f V(/3 2 G).VG < C [ 2 \G\ 2 \F\. 

JDxT 2 JD 

(In this proof G will denote an unspecified constant, changing from line to line.) 
We have 

V(/3 2 G).VG = |V(/3G)| 2 - \G\ 2 \V0\ 2 

so we get 

/ |V(/3G)| 2 <G / |G| 2 VV/3| 2 +/? 2 |G| 2 |F|. (64) 

JDxT 2 JD 

Applying the Sobolcv inequality to \[3G\ and Cauchy-Schwartz we obtain 

||/?G|| 2 6 <G(||G|| 2 2 + ||/3G|| 2 4 ||/|| L2 ). 

(Here we are using the fact that, for fixed K, the metric g is uniformly equivalent 
to the Euclidean metric, so we can transfer the standard Sobolev inequalities to 
our setting.) Using 

\\0G\\U<\\0G\\ L 40G\\l 6 , 

we get 

\\/3G\\h < C (\\G\\ 2 L2 + \\PG\\% 2 \\F\\ L2 \\PG\\% 2 ) <C(£ 2 + £ 3 / 2 ||/?G|| 3 / 6 2 ) . 
This yields 

\\pG\\ L e<C(£ + £ 3 ). 
Thus we have an L 6 bound on G in a neighbourhood of and this gives an 

£3/2 

bound on F * G over this neighbourhood. Going back to (64), we also have an 
L 2 bound on the derivative of (3G: 

l|V(/3G)|| L2 <C(£ + £ 3 ). 

We introduce another cut-off function 7, supported on the neighbourhood where 
0=1. Then we have 

A( 7 |G|) > |V*V( 7 G)| > |A 7 ||G| + |V 7 ||VG| + Cj\F\\G\. 
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Now 

A7 = u lj jij +v l ji. 
So IA7I satisfies a bound, depending on sup \ v\ and we have 

A(7|G|)>*, 

say where 

\\a\\ L 3/2 < C(£ + £ 3 ). 

Then we can apply Theorem 8.15 of [7] (proved by the Moser iteration technique) 
to obtain the desired bound on \G\ at 0. 

For the second part of Proposition 6 we operate with the equation (62). 
We take the L 2 inner product with j 2 F and integrate by parts. For the term 
involving the second derivatives of G we write 

/ 7 2 F.V"V'G= / (V")*( 7 2 F).V'G. 

JDxT 2 JDxT 2 

This yields 

/ |V( 7 F)| 2 <c/ 7 2 |F| 3 + |^| 2 |V 7 | 2 +7|V(7F)||VG + 7|V7||^||VG|. 

JDxT 2 JDxT 2 

Then, using the Sobolev inequlality as before and re-arranging we get an L 2 
bound on the derivative of 7 F near which, in dimension 2, gives the required 
L p bound (since we can suppose that 7 = 1 on the disc \D). 



5.5 Rescaling sections 

In this subsection we will bring together the three ingredients established above 
to obtain a pointwise bound on the tensor G over compact subsets of Q,. This 
involves rescaling the geometric data. A rescaling argument of this kind, using 
balls determined by the Riemannian distance function, would be fairly standard. 
However there are difficulties in carrying this through unless one can establish 
some control of the injectivity radius, or something similar. We get around 
this difficulty by using Caffarelli's theory of the "sections" of a convex function, 
these taking the place of geodesic balls. 

Recall from (5.1) that if u is a smooth, strictly convex function on an open 
set U C R" and x,y is a point in U we have defined H x (y) > 0, vanishing if 
and only if x = y. For t > the section S x (t) at x and level t is the set 

S x (t) = {yeU:H x (y)<t}. 

We will use three results about these sections, or equivalently the functions H x , 
taken from [9] . For each of these results we suppose that the determinant of the 
Hessian satisfies upper and lower bounds 

< A < det(uy ) < A, 
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throughout U, and the constants a below depend only on A, A. Recall that a 
convex set K in R™ is normalised if 

a n B n c K c B n (65) 

where a n = ra~ 3//2 . Any compact convex set can be mapped to a normalised set 
by an affine-linear transformation ([9], Theorem 1.8.2). 

Proposition 7 1. ([9] page 50, Corollary 3.2.4) Suppose that S x (t) is com- 
pact. Then 

Cl t n / 2 < Vol (S x (t)) < c 2 t n/2 . (66) 

2. ([9] page 55, Corollary 3.3.6 (i)) Suppose that S x (t) is compact and nor- 
malised. Then for y € S x (t/2) 

d(y,dS x (t))>c 3 >0. (67) 

3. ([9] page 55, Theorem 3.3.7) Suppose that S x (2t) is compact. Then if 
H x {y) < t and H x (z) < t we have H y (z) < at. 

4- ([9] page 57, Theorem 3.3.8) Suppose that S x (t) is compact and normalised. 
Then S x (t/2) contains the Euclidean ball of radius c$ centred on x. 

The third result can be seen as a substitute for the triangle inequality, if one 
views H x (y) as a defining a notion of "distance" in U. (The assumption that 
S x (2t) be compact does not appear explicitly in [9], where it is assumed that 
U = R™ and all sections are compact, but a review of the proof shows that this 
is the hypothesis needed for our situation.) 

We will now discuss the scaling behaviour of the tensors we have associated 
to a convex function u. For t > and T = Tj e SL(n, R) we set 

u(x) = t^uiVtx) , u*(x) = u{Tx) = t- x u(^tT^x r ) 
We write S = Si for the inverse matrix of T = Tf. Then we have 



det(uy) = det(iiij) = det(u*^); (68) 

\F*\ = \F\=t\F\ (69) 

\G*\ = \G\=t\G\ (70) 

(v*Y - Siv a = VtS{v a . (71) 



Here we have an obvious notation, in which for example v and v* refer to the 
vector fields obtained from the convex functions u,u*. The verification of all 
these identities is completely elementary. 
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With these prclminaries in place we can proceed to our main argument. 
From now on we fix n = 2 and suppose that u is a convex function on ft C R 2 
which satisfies Abreu's equation S(u) = A with constant A. We set 



E= I |F| 2 
In 



p = sup \V\euc- 



Let K CC K + be compact subsets of and suppose that 

< A < det(uy) < A, 
on K + . For x in K we recall that 

H x (dK+)= min H x (y), 



and we let 



We also put 



y£dK+ 



5 = H(K,K+) = minH x (dK+). 

xEK 



D = max H x (dK) . 



Finally, we define a function $ on AT by 

$(a0 = \G(x)\H x (dK), 

and let 

M = max$(a;). 

Theorem 7 There is a constant \i depending only on \,A,E,p,5,D such that 
M <fx. 

To prove this we may obviously suppose that M > 2. We consider a point 
xo where $ attains its maximum value M and set t = |G(xo)| _1 . Then M > 2 
implies that 2t < H Xo (dK) so the section E = S Xo (t) lies in K. Hence S is 
compact in £1 Our first goal is to show that \G\ over S is controlled by i -1 , i.e. 
by its value at xo 

Now we may suppose that 

M > max(2ci4c 4 L'^ 1 ), 

where C4 > 1 is the constant of (3) in Proposition 7, depending on the given 
bounds A, A. This means that if we define e to be 

e = min(l/2,g) 
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we have M > c 2 /e. 

We will now make two applications of the inequality in (3) of Proposition 7. 
For the first we simply observe that in fact S Xo (2t) lies in K and a fortiori in 
K + . Hence this section is compact. Then for y e S Xo {t) we have H Xo (y) < t 
and trivially H Xg (x n ) = < t, so we deduce that H y (x$) < c 4 t. 

Now set r = H X0 (dK) so t = r/M. We claim that if y is in S xo (t) we have 

H y (dK) > -. (72) 

C4 

To see this, the result above yields 

Hy(xo) < cut — C4r/M, 

and c^r/M < re/c 4 since we have arranged that M > c\/e. Set r = re/c^. For 
y € S x „ (t) c K we have 

H y (dK+) >S>^>^ = 2t. (73) 

Suppose that (72) is not true, so there is a point z <E dK with H y (z) < t. Then, 
by (73), Sy{2r) lies in K + and H y (xo) < r so (3) in Proposition 7 would give 

^ ( z ) < C 4T = re < r/2, 

a contradiction to H Xo (dK) = r. 

Now from H y (dK) < re/c^ and the definition of M we obtain our first goal: 

\G(y)\ < C AG(x )\ (74) 

for ally e E = S Xo (t). 

We now invoke our scaling construction for the restriction of u to S. We fix 
the real parameter t to be as above. We know that there is some k > and an 
unimodular affine transformation T so that fci _1//2 T _1 (S) is normalised. But 
we know from (1) of Proposition 7 that the volume of £ lies between c\t and 
c 2 t, so the volume of fct _1 / 2 T _1 (S) lies between k 2 Ci and k 2 c 2 hence 

^ < fc 2 < -. (75) 
8c 2 " " ci v ; 

Thus the convex set S* = t _1/,2 T _1 (S) differs from its normalisation by a scale 
factor which is bounded above and below, so we can apply the results of (2) and 
(4) in Proposition 7, with a change in the constants depending on the above 
bounds for k (alternatively, we could normalise S* by changing the definition 
of t). 
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Lemma 8 There is a constant C depending only on A, A such that 



\v*\euc < C|w| £uc Diam(S). 

To sec this we may suppose that T is diagonal, with = A, T| = A -1 , T 2 = 
T\ = 0. Then, by (71), v* has components 

(v*) 1 = VtA-V , (v*) 2 = VtXv 2 . 

But E* contains the disc of radius of radius R = k~ 1 2~ i / 2 so E = v^^XE*) 
contains an ellipse with semi-axes ^/tXR, VtRX^ 1 . In particular this ellipse is 
contained in £1 so 

2R 

Thus 

K| EUC < Heuc . (76) 

The scaling behaviour (70) of the tensor G and the result (72) established 
above means that \G*\ < c 4 /a over E*. We now apply Proposition 5 to an 
interior set E|5 C E*. There is no loss of generality in supposing that x = 
and that u vanishes to first order at this point, so E* = {y : u*(y) < 1}. We 
define 

SS = {y:u*(y)<3/4}, 

and we set Uq = u* — 3/4. The lower bound on the distance to 9S* furnished by 
(2) of Corollary 7 gives, in an elementary way, a bound on the derivative of Ug 
over Eq and we can apply Proposition 5 to tij, taking D — Eq. Since we have 
bounds on \G* \ and |w*|euc (the latter by Lemma 8), we get an upper bound on 
the Hessian of u'q over the further interior set 

Sl 1 = {y:«*(y)<l/2}. 

The bound on det(u*j) then yields a lower bound on the Hessian over E^. 
We then use (4) of Proposition 7 to see that E!Li contains a disc A about of 
fixed radius, and we can apply Proposition 6 to conclude that 

|G*(0)| 2 < C(£ + £ 3 ), 

where 



/ |F*| 2 < C [ \F* 



(Here we use again the bound on |w*|e uc -) Now the scaling behaviour of F 
implies that 



f \F*\ 2 =t f \F\ 2 <tE. 
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So we conclude that 

|G*(0)| < C(VtE + t 3 E 3 , 

for some constant C depending only on A, A, 5, D, p. But by construction \G* (0) | = 
1, so we obtain a lower bound on t, thus an upper bound on |G(xo)| and thence 
on M. This completes the proof of Theorem 7. 

Given Theorem 7 it is straightforward to complete the proof of Theorem 1 
in the case when A is a constant and Q, C R 2 is a polygon. By the discussion 
of uniform convexity in (5.1) above we can, for any compact set K C fl, find 
further compact sets 

K CC K CC K+ 

in such that 5 = H(K, K + ) and Sq = H(Kq, K) arc bounded below by positive 
bounds depending continuously on the data. Then Theorem 7 implies that the 
tensor G is bounded on K. Covering K by a finite number of sections and 
applying the same argument as above we get upper and lower bounds on the 
Hessian (v,ij) over a neighbourhood of K. Further, we can apply the second 
part of Proposition 6 to get bounds on the L p norm of the tensor F for any 
p. Since tiy is bounded above and below these are equivalent to LP bounds 
on the second derivatives of the matrix (u lJ ). In dimension 2, L\ functions are 
continuous, so we deduce L\ bounds on the inverse matrix (tty). Thus we have 
uniform L\ — hence C 3 ' a — bounds on u over a neighbourhood of K. From this 
point elementary methods suffice to bound all higher derivatives. 
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